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Abstract 

In this paper we dynamically determine the quadrupole mass mo- 
ment Q of the magnetic white dwarf WD 0137-349 by looking for de- 
viations from the third Kepler law induced by Q in the orbital period 
of the recently discovered brown dwarf moving around it in a close 
C^ - ■ 2-hr orbit. It turns out that a purely Newtonian model for the orbit 

of WD 0137-349B, assumed circular and equatorial, is adequate, given 
the present-day accuracy in knowing the orbital parameters of such a 
binary system. Our result is Q = (— 1.5 ± 0.9) x 10 47 kg m 2 for i = 35 
deg. It is able to accommodate the 3-sigma significant discrepancy of 
(1.0±0.3) x 10~ 8 s -2 between the inverse square of the phenomenologi- 
■ cally determined orbital period and the inverse square of the calculated 

Keplerian one. The impact of i, for which an interval Az of possible 



O 



values close to 35 deg is considered, is investigated as well. 
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X ' 1 Introduction 



The recently discovered binary system WD 0137-349 [17, 4, 14], composed 
by a brown dwarf of 0.053 solar masses orbiting a white dwarf of 0.39 solar 
masses along a ~ 2-hr circular orbit, offers us a nice opportunity to dy- 
namically determine the quadrupole mass moment Q of the white dwarf by 
suitably analyzing the orbital period of their relative motion. 

Theoretical calculation of quantities related in various ways to such a 
bulk property of white dwarfs, whose knowledge is important for the equa- 
tion of state of matter at very high densities as those present in such compact 
objects, can be found, e.g., in [6, 7, 16, 1, 11, 8]. 

The approach followed here has been recently adopted to determine Q 
in the double pulsar PSR J0737-3039A/B system [9] and to put upper limits 
on it in binary systems hosting millisecond pulsars [10]. 
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Table 1: Relevant orbital parameters of the WD 0137-349 system [14]. The 
mass of the white dwarf is mi, the mass of the brown dwarf is m 2 and 
the baricentric semimajor axis is denoted with a\ iC . The inclination i is the 
angle between the plane of the sky, perpendicular to the line-of-sight, and 
the orbital plane. 
P h (days) a bc sinz (R Q ) m 1 (M Q ) m 2 (M©) 

0.0803 ± 0.0002 0.375 ±0.014 0.39 ± 0.035 0.053 ± 0.006 



2 Model of the orbital period of WD 0137-349B 

Given the distance and mass scales involved in our problem (see Table 1 for 
the relevant orbital parameters), the first post-Newtonian approximation 
is quite adequate to describe the motion of an object like WD 0137-349B 
around its parent white dwarf. The equation of motion, in standard post- 
Newtonian coordinates, can be written as [23] 

^ = -VU + 1 [-{(3 + 7 ) W 2 + 2( 7 + l)v(v ■ V)U - -fv 2 VU] . (1) 

We will assume that WD 0137-349 rigidly rotates and is endowed with 
both axial symmetry about the equator, taken to be the reference {xy} 
plane 1 . Thus, the gravitational potential U can be written as 



with [18, 12, 22] 



U = U + U Q , 



U = 



tj — _GQ ( 3cos 2 6>-l A 
Q "r 3 " \ 2 J 



(2) 



(3) 



In eq. (3) M = mi + m 2 and 6 is the co-latitude angle (6 = n/2 for points 
in the equatorial plane). The relative acceleration due to the gravitational 
potential of eq. (3) is, in spherical coordinates 





_ _ GM 
r' 2 


< Ag 






r i 




= 0. 



2 r 4 



(4) 



1 Accounting for deviations from such approximations is not needed, given the present- 
day modest accuracy in determining Q. 
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We will now make the simplifying assumptions that the motion of WD 
0137-349B occurs in a circular orbit of radius vq in the equatorial plane of 
the white dwarf. This hypothesis is quite reasonable because tidal forces are 
very strong in short period binaries, and act to quickly circularize the orbit. 
In this case, Ag = = and only the equation for the radial acceleration 
survives in eq. (4) as 

. 2 GM 3 GQ 
r W =^r-- 2 ~ZA- ( 5 ) 

'0 z '0 

By posing ro = a, a simple integration of eq. (5) yields 

Note that eq. (6) is an exact result. 

In regard to the post-Newtonian term p( PN ) coming from the c -2 part 
of eq. (1), in general relativity {(3 = 7 = 1) it is 

p(PN) _ pCl/c 2 ) +p (Q/c 2 ) ) ^ 
where 

pt 1 ^ 2 ) was calculated by Soffel in [20] and Mashhoon et al. in [13]; P^/ c ^ 
can be worked out from (6a) of [19] in the case of equatorial and circular 
orbits. 

Let us now check if the precision with which the orbital period of WD 
0137-349B is known requires to account for the post-Newtonian terms as 
well. Prom Table 1 and by using the following relation for the relative 
semimajor axis a 

1 + — )a hc + 0(c- 2 ), (9) 
mi J 

it turns out 

( P(V C 2 ) = 2 x 10- 7 d, 

(10) 

p(Q/c 2 ) = _ L g x 10 -54q d _ 

Since the uncertainty in the brown dwarf's orbital period amounts to 2x 10 -4 
d, it is clear that the Newtonian model of eq. (5) is quite adequate for our 
purposes. 
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3 Determination of Q 



In order to calculate Q and assess in a realistic and conservative way its 
uncertainty, let us make the following considerations. We are looking for a 
deviation from the 'pure' third Kepler law induced by Q, so that the values 
of the system's parameters to be used with eq. (6) should not come from 
the third Kepler law itself. This is just the case, apart from the inclination 
angle i. Indeed, the mass of the white dwarf comes from an analysis of 
its optical spectrum, independently of any dynamical effect involving the 
orbital motion of its companion [2, 5, 14]; the same also holds for the brown 
dwarf's mass which is determined from the ratio of the masses measured 
by means of accurate spectroscopy of the H a emission and absorption lines 
of WD 0137-349 [14]. Moreover, the dynamical observables at our disposal 
are the orbital period and two semi-amplitude velocities from which the 
projected semimajor axis is phenomenologically determined. The inclination 
i, instead, can only be measured from the expression of the mass function 
T obtained with the third Kepler law, i.e. 

j- _ ™l sin 3 i 
(mi+m 2 ) 2 ' 

according to [14], i 35 deg. Thus, we will not treat i as an estimated 
parameter with an associated experimental error; instead, we will keep it 
fixed to given values close to 35 deg, and for them we will determine Q, 
so to have a realistic idea of what could be the impact of i on Q. We will 
assume Ai/i = 28%, i.e. i = 35 ± 5 deg, and Ai/i = 11%, i.e. i = 35 ± 2 
deg. 

Let us start with i = 35 deg. The values of Table 1 and eq. (6) yield 

Q = -1.4615 x 10 47 kg m 2 . (12) 

According to Table 1 and eq. (6), the error can be conservatively evalu- 
ated as 

SQ < 5Q\ a + 5Q\ M + 5Q\ Ph + 6Q\ G = 0.9004 x 10 47 kg m 2 , (13) 
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with 



SQ\ a < 



2„2| 



5a = 7.401 x 10 46 kg m 2 , 



5Q\m <\ia 2 \5M = 1.449 x 10 46 kg m 



6Q\ Ph < 



SQ\g < 



16 _2 a J 
3 GR? 



2 a 5 



(5P b = 1.49 x 10 45 kg m 2 . 



<5G = 4 x 10 43 kg m z , 



We used SG = 0.0010 x lO" 11 kg~ 4 m 3 s~ 2 [15] and 

5a < 5a\ abc + 6a\ mi + 5a\ m2 = O.O45R , 

with 



6a\ abc < 



Sa\m 2 < 



<Ja| mi < 



1+ m2 
1 mi 



5a hc = O.O27R , 



«]>, 



mi 



(5m 2 = O.O1OR , 



n — 



<J mi = O.OO8R . 



Now, let us repeat the same process for i = 30 deg; we obtain 
Q = (-3.9541 ± 1.8210) x 10 47 kg m 2 . 
For i = 40 deg we get 

Q = (-4.657 ± 4.983) x 10 46 kg m 2 . 



(14) 



(15) 



(16) 



(17) 



(18) 



Let us now see what happens for values closer to i = 35 deg. A departure 
of 2 deg yields 



Q = (-2.1852 ± 1.1754) x 10 47 kg m 



for % = 33 deg and 



Q = (-9.582 ± 7.022) x 10 46 kg m : 



(19) 



(20) 



for i = 37 deg. Thus, for Ai/i = 28% we have AQ/Q w 200%, while a 
narrower variation Ai/i = 11% yields AQ/Q w 83%. Such results have not 
to be considered as experimental errors. 
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4 Discussion 



In regard to the adopted method, let us note the following remarks. 

• The orbital period Pb was determined in a purely phenomenological 
way from spectroscopical measurements, independently of any gravi- 
tational theory, so that it fully accounts for all the dynamical features 
of WD 0137-349B's motion, within the errors. 

• There is a significant discrepancy 

A = (1.0 ± 0.3) x 10~ 8 s~ 2 (i = 35 deg) (21) 

between the inverse square of the phenomenologically determined or- 
bital period and the inverse square of the purely Keplerian period 
1/P( ) 2 which is not compatible with at 3-sigma level; now, eq. (12) 
and eq. (13) yield from eq. (6) 

-|^r = (1.0 ± 0.3) x 10~ 8 s~ 2 (i = 35 deg). (22) 

The same holds also for i = 30 deg and i = 40 deg. 

• According to eq. (12), 3Q/2Ma 2 = -1.2 (i = 35 deg), so that the use 
of the exact, non-approximated expression of eq. (6) is fully justified 

• Let us note that with the result of eq. (12) the second post-Newtonian 
term of eq. (8) becomes p(<2/ c ) = 3 x 10~ 7 d (i = 35 deg): thus we 
can well justify, a posteriori, our choice of neglecting it in our analysis. 

In conclusion, we have determined the quadrupole mass moment of the 
system WD 0137-349 to be Q = (-1.5 ± 0.9) x 10 47 kg m 2 for i = 35 deg. 

Our measured Q should, in fact, be regarded as an effective quadrupole 
mass moment which may also include, in principle, contributions from the 
other multipole mass moments of higher degrees, from the white dwarf's 
magnetic field [11, 8] and from the oblateness of the brown dwarf itself. 

Another possible approach to the problem tackled here would be to re- 
analyze the raw data of the WD 0137-349 system by fitting them with a new 
orbital model including a quadrupole mass term as well, but it is beyond 
the scope of the present work. 

Finally, in order to make easier a comparison with our results, in Table 
2 we quote the numerical values used for the relevant constants entering the 
calculation. 
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Table 2: Values used for the denning, primary and derived constants 
(http: / / ssd.jpl.nasa.gOv/7constants#ref ) . 



constant 


numerical value 


units 


reference 


c 


299792458 


m s _1 


[15] 


GMq 


1.32712440018 x 10 20 


m 3 s~ 2 


[21] 


G 


(6.6742 ±0.0010) x 10" 11 
6.95508 x 10 8 


kg -1 m 3 s~ 2 


[15] 


R© 


m 


[3] 


1 mean sidereal day 


86164.09054 


s 


[21] 
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